Non-extensive random walks 
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The stochastic properties of variables whose addition leads to g-Gaussian distributions Gq{x) 



211/(1-9) 



(with g £ R and where [/(x)]^ = max{f{x),Q}) as limit law for a large num- 



ber of terms are investigated. These distributions have special relevance within the framework of 
non-extensive statistical mechanics, a generalization of the standard Boltzmann-Gibbs formalism, 
introduced by Tsallis over one decade ago. Therefore, the present findings may have important 
consequences for a deeper understanding of the domain of applicability of such generalization. Ba- 
sically, it is shown that the random walk sequences, that are relevant to this problem, possess a 
simple additive-multiplicative structure commonly found in many contexts, thus justifying the ubiq- 
uity of those distributions. Furthermore, a connection is established between such sequences and 
the nonlinear diffusion equation dtp = d^^P" [y 7^ !)■ 



Throughout the last decade increasing attention is be- 
ing given to "non-extensive statistical mechanics" (NSM) 
[1], a theory that extends the standard Boltzmann- 
Gibbs one as an attempt to embrace meta-equilibrium 
or out-of-equilibrium regimes. The cornerstone of the 
non-extensive formalism is the entropy Sq ~ k{l — 
J dx[p{x)]'^)/{q — 1) [2], where fc is a positive constant 
and p a normalized probability density function (PDF) 
(the standard entropy is recovered in the limit q ^ 1). 
Optimization of Sq, under natural constraints, leads to 
q-Gaussian s PDFs defined as 



Gq{x)^Gqmi + {q-l)x^ 



(1) 



where g g M, Gg(0) is a normalization constant, and 
the subindex -f indicates that Gq{x) = if the expres- 
sion between brackets is non-positive. The usual Gaus- 
sian distribution is recovered in the q ^ 1 limit. The 
PDF defined by Eq. (1) also contains the Cauchy and the 
Student's-t distributions, for q~2 and q — (n+3)/(ri+l) 
(where n e N is the number of degrees of freedom), re- 
spectively. 

The extensive research about this proposal brought to 
the surface the fact that many phenomena can be well 
described by g-Gaussian PDFs, from Bose-Einstein con- 
densates [3] or defect turbulence [4] in physics to stock 
returns [5] in finances (see [1] for other examples). De- 
spite the apparent ubiquity of g-Gaussian PDFs, giving 
indirect support to the applicability of NSM, a compre- 
hensive description in terms of first principles is still at 
work. Nevertheless, the fact that g-Gaussian PDFs are so 
frequently found leads to think that a statistical mech- 
anism, as for other limit distributions, may be behind. 
Then, assuming that sums of certain random variables 
fall within the basin of attraction of a given g-Gaussian 
PDF, a basic question arises: Which is the nature of such 
random variables? An answer to this question may give 



useful insights on the general framework of applicability 
of NSM. 

Since sums of independent (or, more generally, weakly 
correlated) random variables lead to either the Gauss or 
Levy limits [6] , depending on the properties of their sec- 
ond moments, then the kind of stochastic variables that 
concern NSM are expected to have some sort of strong 
dependence. This idea is consistent with the nature of 
the systems for which NSM is supposed to apply, e.g., sys- 
tems with long-range interactions or long-term memory. 
To visualize the particular kind of dependence involved, 
it may be instructive to think in terms of one-dimensional 
random walks, since the position of a walker at each in- 
stant is given by the summation of random variables, 
the individual step lengths. As a counterpart, one can 
also investigate the associated evolution equation for the 
probability density of the position of walker. 

It is known that the nonlinear differential equation 

dtp = ^dl.p-- , (2) 

with v, D G M, has as time-dependent solution, for 
p{x, 0) = d{x) and natural boundary conditions [7], 

py{x,t) = G2-u{x/ay{t)), (3) 

with G2-^(0) = g^{t) = l/[7^a^(i)], a^{t) = + 
Dt]-^ and 



r r(^) 




(4) 



r(^) 



The PDF defined by Eq. (3) is normalizable for v > —1 
and has finite second moment for v > 1/3. For v > (<)0, 
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it must be D > (<)0 (see [8] for details). The Cauchy 
distribution corresponds to = 0, in which case there 
is no diffusion. However, in the hniit z/ ^ 0, diffusion 
can be restored by making D oo in such a way that 
Dv > remains finite. The scahng of x with time in Eq. 
(3) indicates superdiffusivc behavior for — 1 < < 1, 
normal diffusion for = 1 and subdiffusion for ;/ > 1, 
although the second moment is finite only for u > 1/3. 

The nonlinear Fokkcr-Planck equation (FPE) (Eq. (2) 
plus a drift term) describes a variety of transport pro- 
cesses, such as percolation of gases through porous media 
[9], thin liquid films spreading imdcr gravity [10] or spa- 
tial diffusion of biological populations [11]. L. Borland 
has shown that it leads to an associated Ito-Langevin 
equation (ILE) with multiplicative noise, where the noise 
amplitude is controlled by a function of the density p [8] . 
In the free-particle, the ILE reads 

x{t) = \D\i[p4x{t),t)]'^r]{t), (5) 

where {rj} is a zero-mean (5-correlated Gaussian process. 
This appears to be a very special kind of multiplicative 
noise that depends on the macroscopic density pi, that, in 
turn, satisfies Eq. (2). Accordingly, the nonlinear FPE 
has been shown to be derivable from a master equation 
where the transition probabilities have the same power 
law dependence on the density p [12]. In terms of tran- 
sition rates, this means that the probability of transition 
from one state to another is determined by the probabil- 
ity of occupation of the first state. This kind of behavior 
is found in many systems ranging from diffusion in disor- 
dered media [9] to financial transactions [5] as discussed 
in Ref. [12]. Notice that, in the limit v ^ 1, one recovers 
the standard Gaussian diffusion, with constant diffusion 
coefficient and transition rates. 

After discretization of time in Eq. (5), the position Xn 
at time t = nr evolves according to 

Xn+i =Xn + v^^[p.(a;„)]('^-i)/2^„, (6) 
where the stochastic process {^}, obtained from ^„ = 
Xi"^^^^ '7(0d<, is Gaussian with (^j) = and {^i^j) = 
5ij . Eq. (5) represents a one- dimensional random walker 
with a sort of statistical feedback in which the steps de- 
pend on the local density of an ensemble of identical 
walkers. Fig. 1 shows the PDFs associated to the ran- 
dom walkers defined by Eq. (6), for different values of 
I/. Notice that for v > 1 a. cut-off occurs at the points 
X = zLa^/y^iy — 1, as a result of the definition of the 
g'-Gaussian in Eq. (1). Convergence is faster as v ap- 
proaches one. 

Let us rewrite Eq. (5) by substitution of the explicit 
expression for pi,{x,t) given by Eq. (3): 

1/2 



X{t) 



1 + (1 - I^) 



al{t) 



riit). (7) 



When V € (—1, 1], the expression between brackets is al- 
ways positive, thus, there is no cut-off. This instance cor- 
responds to long-tailed PDFs, the case of interest in most 



applications. In this regime, the multiplicative noise can 
be interpreted as a unique effective noise for the linear 
stochastic equation 



x{t) = A,{t)r)A{t) + x{t)M,{t)r)M{t), 



where 



M^t) = ,/{l-iy)\D\[g{t)]'^/a4t) , 



(8) 



(9) 



and rjA, tJm arc two independent zero- mean 5- 
correlated Gaussian noise sources with {r]A{t)r]A{t')) = 
{VM{t)VM{t'))=S{t-t') [13,14]. 
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FIG. 1. Distribution of probabilities for the position of ran- 
dom walkers defined by Eq. (6) together with Eq. (3), at two 
times indicated on the figure, (a) v = 4.0 (subdiffusion), (b) 
V — 1 (normal diffusion), (c) v = 0.5 (supcrdiffusion) and (d) 
v — —0.1 (hyperdiffusion). In all cases \D\ = 1 and r — 10~^. 
Symbols correspond to the histogram built from 10^ random 
walkers, starting at xo = (with the addition of a small noise) 
and solid lines to the analytical expression given by Eq. (3). 
Curves collapse after scaling x by t^^^'^'^^K Insets: log-log 
representation to watch the tails. 
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In fact, both Eqs. (7) and (8) yield the same non-hnear 
diffusion equation (2) [14]. Therefore, we have found an 
alternative interpretation for the nonhnear FPE and its 
associated ILE. Notice that, instead of two uncorrelated 
noises, one might consider a unique source with an ap- 
propriate time lag. 

The discretized version of Eq. (8) reads 



2^71+1 — ■^v,n£,A,n + (1 + Mi,^n£^M,n)Xn 



(10) 



where 



M,,„ = ^(l-u)/[\u\iu + l)n\ = sfiilvh , (11) 

and the independent processes {S,a} and {Cm} are zero- 
mean 5-correlated Gaussian ones such that {^A,i^A,j) = 
{^M,iS.M,j) = Sij; as comes out after taking ^x,n = 
Vxmt and = X,(nr)V^, with X G 

{A,M}. The rightmost identities in Eq. (11) define A,y 
and /i, respectively. 

For — 1 < 1/ < 1, the walkers given by Eqs. (6) and (10) 
are equivalent. In fact, it was numerically verified that 
Eq. (10) yields the same PDFs shown in Fig. l(c)-(d), 
which were built from Eq. (6). These PDFs are stable 

through scaling by n'^ . Notice that while A,^^n increases 
with n, My,„ decreases. In the particular case ly = I, one 
recovers the purely additive random process, with a time- 
independent coefficient A,^, i.e., Xn+i = Xn + A^^n- 

The linear charact(T of expression (10) provider the ad- 
vantage of obtaining, by recurrence, an explicit formula 
for the position of the walker, in terms of the random 
sequences {^a} and {^m}, namely, 

n 

Xn+1 = Xo JJ(1 -I- AU^i^Md) + 
n n 

n i^+^-AM,i)- (12) 

In the case v > 1, nonlinearity subsists due to the cut-off 
condition and, therefore, it is not possible to obtain a 
closed form expression as above. When = 1, the prod- 
uct equals one (since for all i), coefficients Ai,j 
are constant, and Xn scaled by n^/^ tends in distribu- 
tion to a stable normal distribution, in agreement with 
the standard central limit theorem. However, as we have 
seen, in the cases ^ I, the distribution of the sum Xn, 
scaled by n^/^"'^^^, tends to a stable (2 — v)- Gaussian. 

The position of the walker at time nr, Xn, can also be 
written as the sum Xn = xq + X^jLo ' where the step 
length at time jr. 



X^ + l - Xj 



Av,i^AA + Mvi£^M,jXj 



(13) 



is somewhat modulated by the full history of the walker. 
Notice that the increments have also a sort of memory 



of the initial condition, through Xj (see Eq. (12)). For 
simplicity, in what follows, we will take xq = 0. Recalling 
that the variables {^a} and {^m} are uncorrelated, and 
noticing that {sj) = 0, then, the two-time correlation of 
the increments is 



SjSji / 



5,r{Al,+Ml^{x])), 



(14) 



where, from Eqs. (11) and (12) 



{x])^Y.'^i,k n i^+Ki) 



fc=0 i=k+l 



r{k + ii + i) 



" (i - D! ^ ■ 



(15) 



By means of the Stirling approximation, for large j 
and v > 1/3, one obtains (x'j) = a1{jT)/{2>u — 1). 
This expression coincides with {x"^) = J x'^Pv{x,t)dx = 
al{t)/{3iy - 1), where t = jr. For ly < 1/3, the Stirling 

formula yields (x|) = (/)(i/)Q!^(_jV)j , where is 

a function of p only, that goes to 1 when -I- 1) ^ 
and diverges at ^/ = 1/3. This result is compatible with 
the fact that (a;^) is divergent in this case, as there is an 
the additional increase of {x"^) / ci^ijT) given by the factor 
with positive power of j. The divergence is logarithmic 
at the marginal case v = 1/3. 

In particular, the increments Sj are ^-correlated, in ac- 
cord with a Hurst exponent H = 1/2 previously found 
[8]. For large j, Eq. (14) leads to 



{iy + l){3,y-l)j 



(16) 



for v > 1/3, and a divergent second moment for u < 1/3. 
It is easy to veriiy, in the former case, that (x^) = 
^J^Q^(s|), as it must be for a sum of uncorrelated vari- 
ables. Moreover, maXj[{Sj)/{xD] ~ 1/n, thus, none of 
the steps gives a predominant contribution to the disper- 
sion of the sum. 

However, higher order moments of the steps do not 
generically factorize. For instance, {s^Sf) 7^ 
Therefore, the steps are not mutually independent al- 
though uncorrelated. This is why the present sums do 
not converge in distribution to the Gauss nor Levy lim- 
iting laws. 

Perhaps, more interesting, are akin processes with 
time-independent coefficients. They are associated to a 
nonlinear FPE with linear drift and yield steady states 
of the g-Gaussian form [14,15]. Namely, the stochastic 
Ito-Langevin equation 



= -7U + ArjAit) + uM rjMit) 
= -ju + [A^ + M^u'^]^r]{t), 



(17) 



with 7 € R, A, M positive constants, and the processes 
{va}, {vm} and {77} defined as above, is associated to 
the following FPE [14] 
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dtp = jdu[up] + ^dlJ[A^ + M^u']p) 



(18) 



Its stationary solution is p1{u) = G2-v{u/ Pu), with 
V = 7/(7 + M2) <l, (3^ = y/T^A/M and G2-^(0) = 
gu = l/ifiulv)- The FPE (18) can be cast in the form 



dtp^^du{up) + %dU[pi{u)r-^p) 



(19) 



where \D\ = A^gl^'^. Then, in this case, the diffusion 
coefficient depends on a power of the long-time solution. 
Time discretization in Eq. (17) leads to 



Un+l = AS^A.n + (1 - K + M^M,n)Un 



(20) 



where the processes {^a} and_{^M} are the same as be- 
fore, K = 7T, A = A^/t and M = My/r. Taking uq = 0, 
and proceeding as for the first walker, the variance is 



A 



n—1 n—1 
k=0 i=k+l 



n A 



^2 ^ ^n-k-l 

k=0 



(21) 



where A — (1— k)^-|-M^. In the limit of large n (and small 
t), (m^) coincides with ^ p\^\u)<lu = 0l/{?>v— 1) — 
AV(27 - M^), for 1/ > 1/3 (hence, A < 1), and diverges 
otherwise (because A > 1). 

Since the stationary distribution is a (2 — j/)-Gaussian, 
this PDF represents the limit law for the random walker 
with constant coefficients. That is, the (2 — i^)-Gaussian is 
the stable distribution of u„ (without further scaling), for 
sufficiently large n. Stability is illustrated in Fig. 1. The 
constant M'^ / k = determines the tail law, through 

the value of while A/M = A/M determines the width 
of the limit distribution. In this case increments 

Tj = Uj+i - Uj = Ai^A,] + (-K + M^M,])uj (22) 

are correlated and history-dependent, as soon as k 7^ 0. 
Because coefficients are not time-dependent, the under- 
lying mechanisms are more evident in this case. 

Summarizing, we have characterized the random vari- 
ables whose addition has a g-Gaussian as limit distribu- 
tion. We considered (i) a purely diffusive situation, where 
the density spreads out, and (ii) a case with external lin- 
ear drift, where a steady density is attained. In both 
cases, the increments originate from two sequences of 
independent, identically distributed Gaussian variables 
{^a} and {£,m}, the first acting additively and the other 
through a linear multiplicative contribution. This struc- 
ture yields a different limit law that the strictly addi- 
tive one that leads to the usual Gauss or Levy lim- 
its. The resulting increments are non- identically dis- 
tributed random variables with null correlation in the 
former case and correlated in the latter. The existence 
of limit laws for similar additive-multiplicative processes 
has been detected before but explicit expressions for the 
PDFs were not found [16]. In conclusion, the ubiquity of 
g-Gaussian PDFs cannot be surprising as far as additive- 
multiplicative structured sequences are in fact frequent 
in several contexts [17]. 
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FIG. 2. Distribution of probabilities for the position of 
random walkers defined by Eq. (20), for three values of 
the number of steps n, indicated on the figure. Parame- 
ters are (a) (A/M, r) = (1,1,0.01), hence f = 0.5; 
(b) {A/M,M'^/j,t) = (1,-11,0.01), hence u = -0.1. Sym- 
bols correspond to histograms built from 10* random walk- 
ers, starting at uo = 0, and solid lines to G2-v{u/ (3^) with 
Pv = \/l — vA/M. Insets: log-log representation to watch 
the tails. 

Following the lines of Ref. [14], the present approach 
could be extended to PDFs of a more general form than 
the g-Gaussian, as for instance p = [1 -f (g— 1)2;^]^^''^ 
with p^2, also frequently found [1]. 

I am grateful to E.M.F. Curado, S.M.D. Queiros, 
A.M.C. de Souza and C. Tsalhs for interesting obser- 
vations. 
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